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Weakly Interacting Massive Particle (WIMP) direct detection experiments are clos- 
ing in on the region of parameter space where neutralinos may constitute the 
Galactic halo dark matter. Numerical simulations and observations of galaxy ha- 
los indicate that the standard Maxwellian halo model is likely to be a poor ap- 
proximation to the dark matter distribution. We examine how halo models with 
triaxiality and/or velocity anisotropy affect exclusion limits, before discussing the 
consequences of the possible survival of small scale clumps. 



1. Milky Way halo modelling 

Current Weakly Interacting Massive Particle (WIMP) direct detection ex- 
periments are just reaching the sensitivity required to detect Galactic dark 
matter in the form of neutralinos. The direct detection event rate and its 
energy distribution are determined in part by the WIMP speed distribu- 
tion, with data analyses nearly always assuming a standard spherical halo 
model with isotropic Maxwellian velocity distribution. 

Observational constraints on the form of dark matter halos depend on 
the relation of luminous tracer populations to the underlying density dis- 
tribution, and are complicated by galactic structure and projection effects, 
however it appears that galaxy halos are triaxial (for a review see Ref. 1 ). 
In particular for the Milky Way, analysis of local stellar kinematics gives 
an estimate for the short-to-long axis ratio of 0.7 ± 0.1 2 , while the great 
circle tidal streams observed from the Sagittarius dwarf galaxy rule out a 
ratio of less than 0.7 in the outer halo 3 (in a flattened potential angular 
momentum is not conserved, so that orbits precess and tidal streams lose 
their coherence). Given the difficulties involved in 'observing' dark matter 
halos it makes sense to turn to numerical simulations for information on 



1 



February 5, 2008 6:9 WSPC/Trim Size: 9in x 6in for Proceedings 



idm02greenrev 



2 



their possible properties. In CDM cosmologies structure forms hierarchi- 
cally; small objects (often known as subhalos) form first, with larger objects 
being formed progressively via mergers and accretion. The shape and in- 
ternal structure of galaxy size halos are then determined by the dynamical 
processes which act on the component subhalos. 

The shape of simulated halos varies, not just between different halos 
of the same mass, but also as a function of radius within a single halo, 
strongly if the halo has undergone a major merger relatively recently. Two 
high resolution Local Group halos studied in detail in Rcf. 4 have axis 
ratios of 1 : 0.78 : 0.48 and 1 : 0.45 : 0.38 at the solar radius. Adding 
dissipative gas to simulations tends to preserve the short-to-long axis ratio 
while increasing the intermediate-to- long axis ratio . 

The anisotropy parameter /3(r), defined as 

< vl > + < v 2 > 

« r >-' \<^> ■ <" 

where < v 2 , >, < v 2 > and < v 2 > are the means of the squares of 
the velocity components, also varies with radius. Typically (3(r) grows, 
although not monotonically, from roughly zero in the centre of the halo 
to close to one at the virial radius, with non-negligible variation between 
halos. The high resolution galactic mass halos studied in Ref. 6 have f3(R@) 
in the the range 0.1-0.4, which corresponds to radially biased orbits. 

2. Effect on exclusion limits 

The differential WIMP event rate due to scalar interactions can be written 
in terms of the WIMP scattering cross section on the proton, a p : 

>o.3 (m p + m x ) 2 2 2 ^ f°° f v 
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where the local WIMP density, p xi is normalised to a fiducial value p .3 = 
0.3GeV/cm 3 such that ( = p x /po.3, A and F(q) are the atomic number 
and form factor of the target nuclei, /„ is the WIMP speed distribution 
in the rest frame of the detector, normalised to unity, and v m - ln is the 
minimum WIMP speed which can produce a recoil of energy E: i> m ; n = 

1/2 

{E(m x + m A ) 2 /2m 2 ,m A ) 

We will consider the two self-consistent triaxial and/or anisotropic 
halo models which have been studied in relation to WIMP direct detec- 
tion to date: the logarithmic ellipsoidal model 7 and the Osipkov-Merritt 
anisotropy model 8 , studied in Ref. 9 . The logarithmic ellipsoidal model 7 
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Figure 1. The speed distributions, in the rest frame of the Sun, for the standard halo 
model (solid line), and the logarithmic ellipsoidal model on the intermediate axis (left 
panel) for parameters which produce axis ratios 1 : 0.78 : 0.48 and (3 = 0.1/0.4 (dot- 
ted/short dashed) and for axis ratios 1 : 0.45 : 0.38 and (3 = 0.1/0.4 (long dashed/dot 
dashed) and for the OM model (right panel) with anisotropy radii which produce 
(3{Rq) = 0.13,0.31 and 0.4 (dotted, short-dashed, and long-dashed). 



is the simplest triaxial generalisation of the isothermal sphere and on either 
the long or the intermediate axis the velocity distribution can be approxi- 
mated by a multi-variate gaussian. Speed distributions on the intermediate 
axis, in the rest frame of the Sun normalised to unity, are plotted in Fig. 1 
along with that for the standard Maxwellian halo model. The logarithmic 
ellipsoidal model has a wider spread in speeds than the standard model, so 
that the differential event rate decreases less rapidly with increasing recoil 
energy. 

In the Osipkov-Merritt (OM) model 8 , which assumes a spherically sym- 
metric density profile, the velocity anisotropy varies as a function of radius 
as j3(r) = r 2 /(r 2 + r 2 ) so that the degree of anisotropy increases with in- 
creasing radius, as is found in numerical simulations. Following Ref. 9 we 
assume a NFW density profile with scale radius r s = 20 kpc. We use values 
of the anisotropy radius r a which correspond to (3(Rq) = 0.14, 0.31 and 
0.4. The resulting speed distributions are plotted in Fig. 1. 

In Fig. 2 we plot the exclusion limits found from the IGEX data 10 , 
taking into account the detector resolution, for the logarithmic ellipsoidal 
model and for the OM model. We also plot the exclusion limits from 
the Hcidclbcrg-Moscow (HM) experiment 11 for the OM model in Fig. 3. 
Comparing Figs. 2 and 3 we see that the change in the exclusion limits 
depends not only on the halo model under consideration, but also on the 
data being used; for IGEX the change in the exclusion limits is largest for 
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Figure 2. The exclusion limits from the IGEX experiment for the logarithmic ellipsoidal 
model, location on the intermediate axis (left panel) and for the OM model (right panel). 



large m x , while for HM the change is largest for small m x . This is because 
for different m x , different energy ranges can be most constraining; for the 
IGEX data the lowest energy bin is always the most constraining, while 
for HM as m x increases the constraint comes from higher energy bins. It 
should therefore be borne in mind when comparing exclusion limits from 
different experiments, that changing the assumed WIMP speed distribution 
will affect the limits from different experiments differently 

The changes in the exclusion limits are not huge (of order tens of per- 
cent) for the experiments we have considered, however these experiments 
are not optimised for WIMP detection. The change in the differential event 
rate, and hence the exclusion limit, would be larger for an experiment with 
better energy resolution and lower threshold energy (see Ref. 12 ). We have 
also seen that different models with the same value for the anisotropy pa- 
rameter [3 can have very different speed distributions, and hence a different 
effect on the exclusion limits. Furthermore it is conceivable that the lo- 
cal WIMP speed distribution may deviate even further from the standard 
Maxwcllian distribution than the models that we have considered. 



3. Small scale clumping 

Numerical simulations arc an extremely powerful tool for understanding 
the large scale structure of the universe, however the local dark matter 
distribution, which is crucial for direct detection experiments, can not be 
probed directly by cosmological simulations. The smallest subhalos resolv- 
able in the highest resolution simulations have mass of order 10 M Q , while 
the first bound neutralino clumps to form would have mass, at matter- 
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Figure 3. The exclusion limits from the HM experiment for the OM model. 

radiation equality, ~ 1O _7 M0 13 . Other approaches 14 A 15 have therefore 
been used to address the problem, with the results obtained depending on 
the method used to extrapolate to small scales below the resolution limit 
of the cosmological simulations a . Moore et al. claim that the phase space 
distribution at the solar radius will depend crucially on the Galaxy's merger 
history and on the internal structure of the smallest subhalos, arguing that 
it is possible that the local dark matter density could be zero, or that a 
single dark matter stream with small velocity dispersion could dominate, 
or that many tidal streams could overlap to give a smooth distribution. 

We will now discuss the consequences for the WIMP direct detection 
rate if small dense clumps survive at the solar radius. We emphasise that 
it is not clear at the moment whether this is the case, and further work is 
necessary to clarify this issue. We could then be currently located within a 
clump with local density in excess of the mean value, on the other hand it 
is possible that we could be in a region between clumps where the WIMP 
density is zero 4 . In the latter case all current attempts at WIMP direct 
detection would be doomed to failure b , and exclusion limits would tell 
us nothing about the WIMP cross-section. At the other extreme a dense 
clump at the earth's location would produce an enhanced event rate and a 
distinctive experimental signal (the theoretical differential event rate would 
be a step function with amplitude inversely proportional to the speed of 
the subhalo with respect to the earth, and position increasing with in- 
creasing relative speed and WIMP mass). For small clump velocities and 



a See Ref. 12 for a more extensive discussion. 

b WIMPs could still be detected using ancient mica though as it has an integtation time 
of order a Gyr 16 . 
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WIMP masses there would be no constraint on the WIMP cross-section (no 
WIMPs would have sufficient energy to cause an observable recoil), but as 
the WIMP mass is increased the constraints would become much tighter as 
then all the WIMPs would be energetic enough to cause events of a given 
recoil energy. 

4. Conclusions 

We have seen that even if the local WIMP distribution is smooth its ve- 
locity distribution may deviate significantly from the standard Maxwellian, 
and this has a non-negligible effect on exclusion limits from WIMP direct 
detection experiments, affecting the limits from different experiments dif- 
ferently. The effect on experiments which can detect the angular and/or 
time variation of the event rate will be more significant. Constraints (and 
in the future possibly best fits) calculated assuming a standard Maxwellian 
halo could be erroneous, even worse if only the signals expected from the 
standard halo model (e.g. a sinusoidal annual modulation with peak in 
June) are searched for, a real WIMP signal could be overlooked. On the 
other hand, more optimistically, if WIMPs were detected it might then be 
possible to derive useful information about the local velocity distribution, 
and hence the formation of the Galactic halo. 
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